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Abstract 

The elastic response is studied of (i) a single flexible chain grafted on a rigid plane and (ii) 
an ensemble of non-interacting tethered chains. It is demonstrated that the entropic theory 
of rubber elasticity leads to conclusions that disagree with experimental data. A modification 
of the conventional approach is proposed, where the end-to-end distribution function (treated 
as the governing parameter) is replaced by the average energy of a chain. It is revealed that 
this refinement ensures an adequate description of the mechanical behavior of flexible chains. 
Results of numerical simulation are compared with observations on uniaxial compression of a 
layer of grafted chains, and an acceptable agreement is shown between the model predictions 
and the experimental data. Based on the analysis of combined compression and shear, a novel 
micro- mechanism is proposed for the reduction of friction of polymer melts at rigid walls. 
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1 Introduction 

This study is concerned with nonlinear elasticity of an individual flexible chain and of an ensemble 
of flexible chains grafted on a rigid surface. This subject has attracted substantial attention in 
the past decade, which may be explained by the importance of the mechanical responses of grafted 
polymer chains and their ensembles in a number of engineering applications, ranging from stick-slip 
transitions in polymer melts near rigid walls at extrusion |^ [2] to drag reduction in dilute polymer 
solutions (reflecting hydrodynamically induced transition from coiled to stretched conformations of 
chains [2|), stabilization of colloidal dispersions by polymers I^, enhancement of adhesion by poly- 
mer brushes and manipulation on polymer membranes and single DNA molecules tethered 
at surfaces [ZllHlini- 

We assume the surface to be sufhciently smooth in the sense that its radius of curvature substan- 
tially exceeds the mean square end-to-end distance of a chain b. Under this hypothesis, a tethered 
macromolecule is thought of as a chain whose end is fixed at some point on a rigid plane, whereas 
all other points are located in a half-space. For definiteness, we suppose that the chain lies in the 
half-space > 0, where a Cartesian coordinate frame {X^} with base vectors {ni = 1, 2, 3) is 
chosen in such a way that the vectors ei and 62 are located in the plane, and 63 is orthogonal to 
the plane. 
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Two models are conventionally employed to describe conformations of a polymer chain. Ac- 
cording to the first, a chain is treated as a random walk with ^ 1 steps that starts at the origin. 
Assuming the length of a step 60 to be constant and small (in the sense that the contour length 
L = boN is finite when N — > 00), we treat the number of steps n as a continuous variable. For a 
flexible chain, the probability density P{n, r) that the nth step ends at a point r is governed by 
the diffusion equation 

where A stands for the Laplace operator, and b = ^Jb^L HTI^. For a Gaussian chain (no restrictions 
on the walk), the function P{n,r) obeys Eq. (1) with the initial condition 

P{0,r) = 6{r), (2) 

where 6 stands for the Dirac delta-function, and the normalization condition 

P(n,r)(ir = l, (3) 



where integration is performed over the entire space. For a grafted chain, the function P{n, r) 
satisfies the additional adsorbing boundary condition jllj 



P(n,r) 



= 0. (4) 

X3=0 



The walk is entirely characterized by the distribution function p(Q) = P{N, Q) of the end-to-end 
vector Q. It follows from Eqs. (1) to (4) that p(Q) is given by 

for a Gaussian chain, and it reads 

for a chain confined to the half-space. Here Qm is the mth Cartesian coordinate of the vector Q, 
Q = IQI, and the subscript index "G" stands for "Gaussian." 

According to the other approach, a chain is treated as a curve (with "length" L) in a three- 
dimensional space. An arbitrary configuration of the chain is determined by the function r(s), 
where r stands for the radius vector, and s E [0, i^]. For a chain that begins at the origin and 
finishes at a point Q, the function r(s) satisfies the boundary conditions 

r(0) = 0, r(L) = Q. (7) 

The energy of a chain is described by a Hamiltonian H{r), which determines the Green function 
(propagator) for the chain 



.r(L)=Q 

G(Q) = / exp 

Jr{0)=0 



H{r{s)) 



V[r{s)], (8) 



where ks is Boltzmann's constant, T is the absolute temperature, and the path integral with the 
measure T>[r] is calculated over all curves r(s) that obey Eq. (7). For a discussion of properties of 
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path integrals, the reader is referred to jl2| As the functional integral is determined up to an 
arbitrary multiplier, an additional constraint is imposed on the function G(Q), 



I G{Q)dQ = 1, (9) 



which allows the Green function G(Q) to be referred to as the distribution function of end-to-end 
vectors p(Q). 

For a Gaussian chain, the Hamiltonian in Eq. (8) reads 

and simple algebra implies that the normalized Green function GciQ) is given by Eq. (5). To 
account for the effect of rigid plane on chain conformations, a penalty functional is inserted into 
the Hamiltonian, 

H{r) = Heir) + ^ ^ 5(r3(s))ds, (11) 



Ljo 

where ^ is a sufficiently large constant, and rm{s) is the mth coordinate of the vector r(s). The 
physical meaning of the integral on the right-hand side of Eq. (11) is the "number" of times that 
the chain crosses the rigid plane in the configuration r(s). When A ^ 1, it is natural to expect that 
the configurations where the chain crosses the plane have a negligible weight in the path integral, 
and they are automatically excluded from the consideration. A rigorous analysis (14) demonstrates 
that this assertion is true: for an arbitrary ^ > 0, the non-normalized Green function reads 

G(».Q) = l-2^e>:p(-^^^),(„)-^ (12) 

with 

ip{x) = —j= — exp(x^)erfc(x), (13) 



ksTh V 2' ' 

where erfc(x) is the complement error function. It is easy to show that after taking the limit A — > oo 
and normalization, the Green function (12) coincides with the function p(Q) given by Eq. (6). 

The classical theory of entropic elasticity has been formulated more than half a century ago, 
and its exposition can be found in a number of monographs, see, e.g., jl5l I16j . This concept is 
grounded on the Boltzmann formula that expresses the probability density p(Q) in terms of the 
free energy of a chain ^'(Q), 

p(Q) = exp(-|M). (14) 

According to the finite elasticity theory, two states of a chain are distinguished: (i) the reference 
(initial) state occupied before application of external loads, and (ii) the actual (deformed) state 
that the chain acquires after deformation. The end-to-end vector of a chain in the deformed state 
Q' is expressed in terms of that in the reference state Q by the formula 

Q' = F-Q, (15) 

where F is the deformation gradient, and the dot stands for inner product (for simplicity, we adopt 
the affinity hypothesis, according to which, the deformation gradient at the micro-level coincides 
with the deformation gradient for macro-deformation). It follows from Eqs. (14) and (15) that the 
increment of free energy A^'(F, Q) driven by deformation of a chain reads 

A^'(F, Q) = -fcBr[lnp(F • Q) - Inp(Q) . 
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The strain energy per chain T4^(F) is calculated by averaging the increment of free energy over the 
initial distribution of end-to-end vectors, 



W{F) = ksT J [inp(Q) - lnp(F • Q)]p(Q)ciQ, (16) 



where the integration is performed over the entire space. Given a strain energy W{F), constitutive 
equations for a chain or an ensemble of chains (whose strain energy equals the sum of strain energies 
of individual chains) can be determined by conventional formulas. For example, for a Gaussian 
chain with the distribution function (5), Eq. (16) implies the classical formula 



W = hBT(li{C)-s), 



where I„i stands for the mth principal invariant of a tensor, C = F""" • F is the right Cauchy-Green 
deformation tensor, and T denotes transpose. For a chain grafted on a rigid plane, the natural 
restriction is imposed on the deformation gradient F: the end-to-end vectors before and after 
deformation, Q and Q', are located in the same half-space. Although this constraint complicates 
calculations in the general case, stress-strain relations can be easily developed in an explicit form 
for deformation programs traditionally studied in nonlinear mechanics. 

The starting point of the present work is the difference between the two models for polymer 
chains. For a chain modeled as a random walk, the Boltzmann equation (14) provides the only 
way to introduce free energy. On the contrary, for a chain described as a curve with a Hamiltonian 
H ascribed to it, the free energy ^'(Q) may be identified as the average value of the Hamiltonian 
{H)q, 

Hiris))- 



^^^^ = rlr^ ,. //(r(.))exp 

<^\^) Jr(0)=O 



kBT 



V[ris)], (17) 



1 rriL)=Q 

G(Q)A(o)= 

where the pre-factor 1/G(Q) plays the role of the partition function in the conventional formula 
for averaging. Given a free energy ^'(Q), we can define the non-normalized distribution function 
of end-to-end vectors p^(Q) by Eq. (14), and, after appropriate normalization of p*(Q) with the 
help of Eq. (3), calculate the strain energy per chain W{F) by the formula similar to Eq. (16), 

W{F) = J [*(F • Q) - *(Q)]p*(Q)dQ. (18) 

Three questions arise naturally after the formulation of this "upside-down" (with respect to the en- 
tropic elasticity theory) approach to the derivation of constitutive equations: (i) does this technique 
lead to classical results for a Gaussian chain, (ii) is there any case of practical interest where the 
novel method implies physically plausible conclusions, whereas the traditional approach fails, and 
(iii) how difficult is it to calculate the path integral in Eq. (17) provided that the Green function 
G(Q) is known? The objective of this paper is to shed some light on these issues. 

The exposition is organized as follows. In Section 2, a method is described that allows the 
average free energy ^(Q) to be found from Eq. (17) without calculation of the path integral. In 
particular, it is shown that for a Gaussian chain, constitutive equations (16) and (18) coincide. 
Force-stretch relations for an individual chain grafted on a rigid surface are developed in Section 
3. It is demonstrated that two approaches lead to qualitatively similar conclusions for tension 
(compression) of a chain, although the results differ quantitatively. Stress-strain relations for an 
ensemble of noninteracting fiexible chains grafted on a rigid surface are developed in Section 4. It is 
revealed that the conventional theory of rubber elasticity fails to adequately describe observations, 
whereas the non-entropic concept provides good agreement with experimental data. Combined 
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shear and compression of a layer of non- interacting grafted chains is studied in Section 5, where a 
novel mechanism is suggested for the reduction of friction of polymer melts near rigid surfaces (in 
addition to the standard entanglement-disentanglement process). Some concluding remarks are 
formulated in Section 6. 



2 Calculation of the free energy 



Our aim is to determine the function ^'(Q) from Eq. (17) explicitly, provided that the Green 
function G(Q) for a flexible chain is known. To simplify transformations, it is convenient to 
consider a more general problem of averaging for an arbitrary Hamiltonian H of the form 



H(y) = Hg(y) + A^{v), 



(19) 



where $(r) is an arbitrary potential that describes intra-chain and inter-chain interactions. Ac- 
cording to Eqs. (8) and (10), the Green function for the Hamiltonian (19) reads 



/•r(L)=Q r i /'ih-oT f^/rlr \2 

OniAMM)-! exp[-^(^/^ (-(.)) d. + ^^(r(. 



V[v{s) 



(20) 



where the parameters 6o and A are included explicitly as arguments of the function Gh- Differen- 
tiating Eq. (20) with respect to 6o and using Eqs. (10) and (17), we find that 



dc 



H 



1 



r(L)=Q 



ksTbo Jr(o)=o 

-{Hg)ciGh- 



i?G(r(s))exp 



i-^H^ / (-ris)) ds + A^(r(s 



P[r(s) 



ksTbo 

Similarly, differentiation of Eq. (20) with respect to A results in 

OGh 1 



(21) 



OA ksT 
It follows from Eqs. (17), (21) and (22) that 



{^)ciGh. 



(22) 



*(Q) 



{Hg + yl$)Q 

keT 
GH{A,bo,Q) 



bo^iAboM) - A—^{A,boM) 



dbo 



dA 



(23) 



Formula (23) provides an analytical expression for the average free energy of a flexible chain with 
an arbitrary potential For example, for a Gaussian chain {A = 0), Eqs. (5) and (23) imply that 



(24) 



Evidently, insertion of expression (24) into Eq. (14) and subsequent normalization results in the 
distribution function given by Eq. (5), which means that for a Gaussian chain, our approach 
coincides with the conventional one. 

For a flexible chain confined to a half-space, we calculate the averages {Hg)ci and {A^)q 
separately and use the first equality in Eq. (23). Equation (21) results in 



{Hg)q 



keTbo dGn 
Gh dbo ■ 
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As the differentiation witli respect to 5o is independent of A, we can use Eq. (6) for the limit of 
the Green function when A — > oo. Substitution of Eq. (6) into this equahty imphes that 



{Hg)q = knT^ 



3Q2 
262 



2 . 



(25) 



Comparison of Eqs. (24) and (25) shows that the average Gaussian Hamiltonian is not affected by 
the presence of the constraint (additive constants in these relations do not influence the increment 
of free energy A^'). 

It follows from Eqs. (12), (13) and (22) that for any ^ > 0, 



Q 



ksTa dG 



(a,Q), 



Gh dA G{a,Cl) da 
where G{a, Q) is given by Eq. (12). Differentiation of Eq. (12) with respect to a results in 



(26) 



dG 
da 



(a,Q) 



exp 



3Qir2 



262(1 -r2 



ip'iar) 



dT 



72' 



where the prime denotes differentiation with respect to the argument. Bearing in mind that 

^(ot) = acp [aT) 
and integrating by parts, we present this equality in the form 



exp 



(- 



3Qir2 



262(1 _r2)yyr—f 
3g2^2 



r=l 



T=0 



d 



^ai ^^""^^d^lM 



262(1 -r2)yyr^ 



dr. 



(27) 



The out-of-integral term vanishes (because of the property of the exponent at r = 1 and due to 
Eq. (13) for the function (p{x) at r = 0). Insertion of expression (27) into Eq. (26) yields 



lim {A^)q = -2kBT 



lim^^oo Jo a¥'(or)|: [exp(-2^^^ 



dT 



lima^oo aG{a, Q) 
The limit of the function aG{a, Q) was calculated in jTl], 



Taking into account that 



lim aG{a,Q) = ^J^. 



lim a(p{aT) 



and performing differentiation with respect to r, we arrive at the formula 



lim (A$)q 



ksTb nr 

Q3 V Sir 

ksTb FT 
Q3 V Sir 



exp 



(- 



3gir2 



f— ^ 
V62(l 



3Q 



dT 



262(1 -r2)yV62(l-r2) 



(1-t2 
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where t = r/Vl — t^- Calculation of the Gaussian integral results in 

Jim(^$)Q = fc.T(^ + l), 
which, together with Eqs. (23) and (25), implies that 

<f(Q) = teT(^ + ||-l). (28) 

Formula (28) provides an explicit expression for the average free energy of a flexible chain grafted 
on a rigid surface. 



3 Uniaxial tension of an individual chain 

Our aim now is to compare the mechanical response of an individual flexible chain grafted on a 
rigid plane when its free energy is (i) determined by the entropic elasticity theory, Eq. (14), and 
(ii) given by Eq. (28). 

We begin with the conventional approach, substitute expression (6) into Eq. (14), and find that 

*.(Q) = feT(g-l„^). (29) 

where the subscript index "e" stands for "entropic." To determine the natural state of the chain, 
that is the state (described by an end-to-end vector Q°) in which the free energy has its extremum, 
we differentiate Eq. (29) with respect to Qm, equate the derivatives to zero, and obtain 

Q? = o, go = o, ^3 = ^- (30) 

It follows from Eqs. (29) and (30) that for a fixed Q3 (shear), the response of a tethered flexible 
chain coincides with that of a Gaussian chain. Thus, we concentrate on tension (compression) 
perpendicular to the rigid plane. Denote by z the displacement (from the natural state) of the free 
end of the chain along the 63 vector when a force f = fe^ is applied to this end. Given an elastic 
energy the force / is expressed in terms of the displacement z by the formula 



(31) 

Q3=Q°3+Z 



dQ 

Combination of Eqs. (29) to (31) implies that 

The first term on the right-hand side of Eq. (32) describes the mechanical response of a Gaussian 
chain, whereas the last term characterizes the influence of the constraint on the force-stretch 
relation. The expression in the parentheses approaches zero as 2 ^ oc, which means that at large 
extensions, the effect of surface becomes insignificant. This expression is important at compression 
{z < 0), because it tends to infinity when z —Q3 (the total compression of the chain requires an 
infinite force). The stiffness (an equivalent spring constant) of the chain fi is determined as 



df 

^ = Tz 



(33) 

2=0 
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Equations (32) and (33) imply that 

/ie = QknTb-^, (34) 

which means that the presence of a rigid surface increases the stiffness of a Gaussian chain hg = 
SkBTb-"^ dZI by twice. 

We now repeat the same calculations for a tethered flexible chain with free energy (28). Differ- 
entiating Eq. (28) with respect to Qm and equating the derivatives to zero, we find the end-to-end 
vector in equilibrium 

g; = o, Q° = 0, Q3 = fe/|- (35) 
Substitution of expressions (28) and (35) into Eq. (31) implies the force-stretch relation 

/■«(^)-^K-^(-(^n/|)l> 

where the subscript index "ne" stands for "non-entropic." Equations (33) and (36) result in 

/ine = UkBTb-^, (37) 

which means that the stiffness of a flexible chain grafted on a rigid surface exceeds that of a Gaussian 
chain by a factor of four. 

Equations (32), (34) and (36), (37) show that, although the force-extension relations derived 
within these approaches are quantitatively different, they are quite similar qualitatively. To reveal 
a qualitative difference between the two concepts, we study the elastic behavior of an ensemble of 
grafted chains. 



4 Uniaxial tension of an ensemble of chains 

We now consider the response of an ensemble of non-interacting flexible chains grafted on a rigid 
plane. In the rubber elasticity theory, inter-chain interactions are conventionally accounted for in 
terms of the incompressibility condition JUj , which means that the neglect of interactions between 
chains is tantamount to the assumption about compressibility of the ensemble. The difference 
between the analysis of an individual chain and that of an ensemble of macromolecules is that we 
do not assume end-to-end vectors of chains in an ensemble to be in their natural states, but suppose 
that the distribution of end-to-end vectors Q is governed by an appropriate probability density. 
At uniaxial tension in the direction orthogonal to the plane, the deformation gradient F reads 

F = eiei 6262 Aeses, (38) 

where A stands for the elongation ratio. It follows from Eqs. (6), (14) and (38) that within the 
entropic elasticity theory, the increment of free energy per chain is given by 

A^e(A, Q) = kBT{^{X^ - 1) - In a). 

Inserting this expression and Eq. (6) into Eq. (16) and calculating the integral, we find the strain 
energy of an individual chain, 

We{X) = kBT{X^ -In X-1). (39) 

Denote by C concentration of grafted chains (the number of chains per unit area of the surface). 
Multiplying the strain energy of a chain W by the number of chains per unit area, we determine 
the strain energy per unit area of a layer, 

W = CW. (40) 
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At uniaxial extension, the work of external forces (per unit area and per unit time) reads 



A dt ' 

where S is the tensile stress, X^^dX/dt is the rate of longitudinal strain, and h is the initial height 
of the layer. According to the first law of thermodynamics. 



dW 
~dt 



n. 



Combining these relations, we arrive at the formula 

S-A^^^ 

^-\l[x- 



(41) 



(42) 



Equation (42) is similar to the well-known formula that expresses the principal Cauchy stresses 
in terms of the principal elongation ratios A^, 



5jrr), — Xr 



dW 



-'m — ^^m J^ ) 



where W is the strain energy per unit volume in the reference state and m = 1,2,3. Inserting 
expression (39) into Eq. (42), we find that 



Se(A) = ^(2A2-l). 



(43) 



Before discussing the physical meaning of Eq. (43), we find a stress-strain relation for a layer of 
grafted chain with the average free energy given by Eq. (28). It follows from Eqs. (28) and (38) 
that the increment of free energy A^'ne reads 



A*„e(A,Q) = ^.r[g(A2-l) + ^(^-l) 



262 ^' ' 3Q2VA2 

The expression for the distribution function of end-to-end vectors is found from Eqs. (14) and (28), 



3Q2 62 



Pne(Q) = C'exp -(jyTT + 



262 3Qi 



(44) 



where the pre-factor C is determined by the normalization condition (3) . Calculation of the integrals 
over Qi and Q2 in this relation implies that 



C 



1 /3 



Cn 



1 



27r-fe)'^°' 

where z = Qs\/3/b. Substitution of these expressions into Eq. (18) results in 

VFne(A) = kBT[Ci{X^ _ 1) + C2(^ - l) 

where 



C, = C,l exp(-(-.2 + _))_,,, c, = C,l exp(-(-.2 + _)) 



1 \\ 

^2 ■ 



(45) 



(46) 



(47) 
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It follows from Eqs. (42) and (46) that 



Sne(A) = 2^(CiA2-C2A-2). (48) 

To find the coefficients Ci and C2 explicitly, we transform the second equality in Eq. (45) by 
integration by parts 

±=2/~exp(-(l.V^))(f-^)<i.. 
Combining this relation with Eq. (47), we conclude that 

Ci - C2 = ^. (49) 



We now set y = ■\/2/z in the second equality in Eq. (45) to obtain 

^ = V2/^~exp(-(i/ + J,))J,„ 
It follows from this formula and Eq. (47) that 



V2 

Equations (48) to (50) yield 

-XkBT 



C2 = (50) 



Sne(A) = V2^[(l + ^)A2-A-^]. (51) 



V2 

According to Eqs. (43) and (51), the tensile stress Sq = CksT/h should be applied to an ensemble 
of grafted chains in order to maintain its initial (undeformed) state. Splitting the longitudinal 
stress S into the sum of the residual stress Sq and the extra-stress a, we find that 

ae(A) = 2^(A2-1), (52) 
MA) = V2^[(l + i=)(A^-l)-(A-2-l)]. (53) 

At strong extension of a layer of chains (A 3> 1), Eqs. (52) and (53) show a monotonic increase in 
tensile stresses with A, but the non-entropic model implies a higher (by a factor of ^(\/2 + 1)) rate 
of growth of the extra-stress a. 

Formulas (52) and (53) demonstrate the principal difference between the two approaches at 
compression. According to the entropic elasticity theory, Eq. (52), a finite compressive stress 
CTg = — 2So leads to the total compression of a layer of flexible chains (A = 0). As this conclusion 
contradicts basic hypotheses of continuum mechanics, the corresponding strain energy density, Eq. 
(39), should be excluded from the consideration. On the contrary, Eq. (53) reveals that an infinite 
compressive stress is required for the total compression of a layer of grafted chains, in accord with 
the axioms of the nonlinear theory of elasticity. 

It is worth recalling that Eqs. (52) and (53) are derived for a compressible layer of grafted 
chains, which implies that they differ from conventional relations describing uniaxial extension of 
an incompressible medium (where the engineering stress is proportional to A, not to A^). This 
discrepancy is driven by the fact that the area of a compressible layer with the normal vector 63 
does not change, while at uniaxial tension of an incompressible material, this area decreases as A~^. 
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It is of interest to compare predictions of the non-entropic model with experimental data. For 
this purpose, we use observations at compression of red cell membranes by a micro-sphere tip (the 
bio-membrane force probe). A red cell ghost after a preliminary treatment was put on a glass 
substrate and was compressed by a glass bead with a micron-size radius. The compressive force 
fc was measured simultaneously with the distance from the bead to the substrate z. The red cell 
membrane on the substrate is modeled as a rigid layer, to which a compressible spectrin network 
(treated as a layer of non-interacting flexible chains) is linked by junctional complexes of short 
actin filaments and other proteins. For a description of the material, the experimental procedure, 
and the method of analysis of measurements, the reader is referred to the original publication 
In numerical simulation, the force fc is found from Eq. (53) where the stress Une is multiplied by 
the contact area S (following P, the latter is assumed to be independent of z). The elongation 
ratio A is connected with the distance z by the formula A = z/Z, where Z = 205 nm is the 
initial distance (determined from the measurements as the smallest distance at which the force 
fc vanishes). The compressive force fc is plotted versus distance from the substrate z in Figure 
1 together with the curve fc{z) calculated at VlCksTS/h = 30 pN. Fi gure 1 demonstrates an 
acceptable agreement between the observations and the results of numerical analysis. We do not 
treat this agreement as a confirmation of the model (the experimental data were obtained based 
on a number of simplifications that may be questioned), but rather as a demonstration of failure 
of the entropic theory of rubber elasticity. 



5 Superposition of uniaxial tension and shear 

Our purpose now is to consider the influence of uniaxial tension (compression) on shear deformation 
of a layer of grafted flexible chains. An exposition of motivation for the study of this problem is 
postponed to the discussion at the end of this section. The analysis is focused on the non-entropic 
model with the average free energy (28), whose results are compared with those of the entropic 
elasticity theory. 

Subsequent imposition of uniaxial tension with an elongation ratio A and the deformation 
gradient Fi given by Eq. (38) and simple shear with the deformation gradient 

F2 = eiei + 6262 + 6363 + ^6263, 

where n stands for shear, induces deformation with the deformation gradient 

F = F2 • Fi = eiei + 6262 + Ae3e3 + AKe2e3. (54) 

It follows from Eqs. (28) and (54) that 

A^nc(A, K, Q) = fcBrj^ [((Q2 + XKQsf - Qi) + (A2 - l)Ql\ + ^(A-' - !)}• (55) 

Inserting expressions (44) and (55) into Eq. (18) and calculating the integral, we find that 

W {X, k) = kBT[Ci[x^ {I + K^)-l)+ 02(^-1)], (56) 

where the coefficients Ci and C2 are given by Eq. (47). At combined uniaxial tension and shear, 
the work of external forces (per unit area of the layer of grafted chains and unit time) reads 
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where Si denotes the (Cauchy) shear stress. The remark that Si is defined per unit area in the 
deformed state is important: at uniaxial tension Si is less than the engineering shear stress by a 
factor of A"-*^, while the shear displacement of the upper boundary of the layer increases by a factor 
of A, which means that no additional multiplier appears in the expression for the work of shear 
forces. Substitution of this expression and Eq. (40) into Eq. (41) results in 



A 



dX 



CdW\ 
h~dx) dt 



+ Si 



CdW\dK 
h dn / dt 



0. 



Bearing in mind that this equality is satisfied for arbitrary functions X{t) and K{t), we conclude 
that the expressions in parentheses vanish. Equating the first expression to zero, we arrive at Eq. 
(42), whereas equating the other expression to zero, we obtain 



Si = ^ 



CdW 



h dn 



(57) 



It follows from Eqs. (42), (56) and (57) that 



h 



Ci(1 + k2)a2-C2A 



-2 



2^Ci.A^ 



Inserting expressions (49) and (50) into these relations and splitting the longitudinal stress S into 
the sum of Sq and the extra-stress cr, we find that 



a = V2' 



h 



1 



V2 



1 + ^(1 + - \- 



1 

71 



(1 + V2)^acA2. 



(58) 



Introducing the notation 



a = 



ah 



r = 



Si^ 



and excluding A^ from Eqs. (58), we obtain 



a 



T- 



1 + 



K 



1 + V2 k _ ^ 
V2 f 



Resolving this equation with respect to r, we find that 

^2 



T = 



2(1 + ^2) 



(l + a)+ (1 + a)2 + 2^2(1 + V2) 



1 + 



(59) 



The dependence of the dimensionless shear stress f on shear k is plotted in Figure 2 at various 
values of the dimensionless tensile stress a. This figure demonstrates that given a, the shear 
stress monotonically increases with shear. Given k, the stress f pronouncedly grows with a. At a 
relatively large compression (a = —5.0 to —20.0), the function f{K) is practically linear, and the 
tangent shear modulus is rather small. When a is positive (tension), the function r(«;) becomes 
strongly nonlinear: it rapidly grows at relatively small k and increases as at large k. 

It is of interest to compare the solution (59) with that found by using the entropic elasticity 
theory. In the latter case, an analog of Eqs. (58) is given by 



a = 2 



h 



(1 + «2^;^2 _ 1 



Si=2^.A^. 
h 



Introducing the dimensionless variables a and r and excluding A^ from these relations, we arrive 
at the formula 

^'^'^^ (60) 



T = 



l + /«2 
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Equation (60) contradicts the physical intuition: it shows that given a tensile stress a, the shear 
stress El depends on shear k non-monotonically: Si increases with /t at small strains, reaches its 
maximum, and decays to zero at large deformations. 

It is necessary to provide some explanations for our choice of shear deformation with a fixed 
normal stress for the analysis. This problem naturally arises in the study of laminar flows of polymer 
melts near a rigid surface. Experiments reveal that some chains from the melt are grafted on the 
surface. At low flow velocities and, as a consequence, small pressures in the melt (in the model, 
the pressure p is equivalent to the compressive stress a), a strong friction is observed between the 
melt and the rigid wall which suppresses slippage of the melt entirely. At higher flow velocities, 
measurements demonstrate a significant reduction in friction and noticeable slippage of the melt 
with respect to the surface QEI; in particular, when the grafted chains do not overlap (an ensemble 
of non- interacting tethered chains). Two reasons for this decrease in friction may be mentioned: 
(i) disentanglement of polymer chains in the bulk from grafted chains ^ |5] , and (ii) a decay in the 
tangent shear modulus of the layer of grafted chains driven by the growth of pressure. Although 
the entanglement-disentanglement mechanism may be dominant, the above analysis demonstrates 
that an increase in pressure p may provide substantial contribution into the decrease of friction as 
well. 

6 Concluding remarks 

The nonlinear elastic behavior has been analyzed of an individual flexible chain and an ensemble of 
non-interacting chains grafted on a rigid surface. It is demonstrated that the conventional entropic 
elasticity theory leads to the conclusions that contradict our physical intuition, which implies 
that its applicability in the mechanics of polymers is questionable. A modification of this theory 
is proposed, where the average energy of a fiexible chain is treated as the governing parameter 
instead of the distribution function of end-to-end vectors. It is shown that this refinement leads to 
physically plausible conclusions for an ensemble of tethered fiexible chains under uniaxial tension 
(compression) and under superposition of uniaxial tension and shear. The results of numerical 
simulation for the former problem demonstrate fair agreement with observations on compression 
of red cell membranes. Our results for the latter problem reveal a novel micro- mechanism for the 
decrease in friction (with the growth of fiow velocity) of polymer melts moving near solid walls. 
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List of figures 

Figure 1: The compressive force fc versus the distance from the substrate z for a red cell mem- 
brane. Circles: experimental data [O]. Solid line: results of numerical simulation. 

Figure 2: The dimensionless shear stress f versus shear k for a = —20.0, —10.0, —5.0, —2.0, —1.0, 
0.0, 2.0 and 5.0, from bottom to top, respectively. 
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